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Although the convergent close-couphng (CCC) method has 
achieved unprecedented success in obtaining accurate theo- 
retical cross sections for electron-atom scattering, it generally 
fails to yield converged energy distributions for ionization. 
Here we report converged energy distributions for ionization 
of H(ls) by numerically integrating Schrodinger's equation 
subject to correct asymptotic boundary conditions for the 
Temkin-Poet model collision problem, which neglects angular 
momentum. Moreover, since the present method is complete, 
we obtained convergence for all transitions in a single calcu- 
lation (excluding the very highest Rydberg transitions, which 
require integrating to infinitely-large distances; these cross 
sections may be accurately obtained from lower-level Rydberg 
cross sections using the 1/n^ scaling law). Complete results, 
accurate to 1%, are presented for impact energies of 54.4 and 
40.8 eV, where CCC results are available for comparison. 



PACS number(s): 
03.65.Nk 



34.80. Dp, 34.80. Bm, 31.15.Fx, 34.10. +x. 



The Temkin-Poet (TP) model ||I|,|| of electron- 
hydrogen scattering is now widely regarded as an ideal 
testing ground for the development of general methods 
intended for the full three-body Coulomb problem. Al- 
though only s-states are included for both projectile and 
atomic electrons, this model problem still contains most 
of the features that make the real scattering problem hard 
to solve. Indeed, even in this simplified model, converged 
energy distributions for ionization can not generally be 
obtained via the close-coupling formalism H. Any gen- 
eral method that can not obtain complete, converged 
results for this model problem will face similar difficul- 
ties when applied to the full electron-hydrogen system. 
Therefore we believe it is essential to develop a numerical 
method capable of solving the TP model completely be- 
fore angular momentum is included. Here we report such 
a method. Complete, precision results for e~ -I-H(ls), ac- 
curate to 1%, are presented for total energies of 3 and 2 
Rydbergs (Ryd). Atomic units (Ryd energy units) are 
used throughout this work unless stated otherwise. 

Our numerical method may be summerized as follows. 
The model Schrodinger equation is integrated outwards 
from the atomic center on a grid of fixed spacing h. The 
number of difference equations is reduced each step out- 
wards using an algorithm due to Poet y, resulting in a 
propagating solution of the partial-differential equation. 
By imposing correct asymptotic boundary conditions on 
this general, propagating solution, the particular solu- 
tion that physically corresponds to scattering is obtained 
along with the scattering amplitudes. 



The Schrodinger equation in the TP model is given by 
2 



\dx'^ dy^ min(a:,y) 



E\^[x,y) = ^, 



with boundary conditions 

*(x,0) = 'J'(0,y) = 
and symmetry condition 

*(y,a;) = ±*(a;,y). 



(1) 



(2) 



(3) 



depending on whether the two electrons form a singlet 
(-I-) or triplet (— ) spin state. Eq. (1^) is separable in the 
two regions x > y and x < y. Because of the symme- 
try condition (|^), we can solve Eq. (|^) in just one of 
these regions and this is sufficient to determine all of the 
scattering information. For brevity, we do not explicitly 
indicate the total spin since the singlet and triplet cases 
require completely separate calculations. For x > y, the 
wave function may be written 



-ik^ . X 



^C,^i^p,^{y)e'"'^"' 



^{x,y) = Vei(2/)e 



The ipe are bound and continuum states of the hydrogen 
atom with zero angular momentum: 



+ / debCe,,V'e.(y)e^'==^^ (4) 



Aiy) = ye"'\F,{l-l/q,2;2qy). 



(5) 



Here q^ = —e, where e is the inner electron energy, and 
li^i is the confluent hypergeometric function. The mo- 
menta in (y) are fixed by energy conservation according 
to 



£b 



E, 
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where £^ > is the total energy. The C^i are related to 
S-matrix elements by normalization factors: 



hi 
for discrete transitions and 
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for ionization, where k = ^Jth- Cross sections are ob- 
tained from S-matrix elements in the usual manner. 

To convert the partial-differential equation (nl) into dif- 
ference equations we impose a grid of fixed spacing h and 
approximate derivatives by finite differences. After ap- 
plying the Numerov scheme in both the x and y direc- 
tions, our difference equations have the form Ej 



Here we have collected the various ^.'^j — f,2,. 



^3 



,(») 



(9) 

where ^j*'' = ^(x — ih, y — jh), into a vector ^ . The 
matrices A^'^ B^'^ and C*^*' are completely determined 
by the formulas given by Poet [El. 

At each value of i we can solve our equations if we 
apply symbolic boundary conditions at i -f 1 [solve for 
\E'^- in terms of \1/^' (j = 1, 2, . . . , i)]. This procedure 
yields a propagation matrix D'*^: 



* 



(■0 



d(*) • * 



■('+1) 



(10) 



We can obtain a recursion relation for D^'^ by using (IT 
to eliminate ^ from equation (^): 

bW + A^'^ ■ D^'-i)] • #'' = -CW • W^*"^^\ (11) 

Comparing ( [ll| ) with (|lO|), 



• c«. 



(12) 



Thus each D^*) is determined from the previous one (D*^^^ 
can be determined by inspection). 

In the asymptotic region, the form of the wave function 
is known and is given in terms of the C^i by 



* 



{i) 



l(*) + r(«) . c. 



(13) 



Here the matrix I*^*' contains the incident part of the 
asymptotic solution while R^*^ contains the reflected 
part. The asymptotic solution is identical to the full 
solution, Eq. (0), except that the quadrature over the 
continuum extends only up to the total energy E. The 
infinite summation over discrete channels is truncated to 
some finite integer Nd and the quadrature over the two- 
electron continuum is performed prior to matching by 
first writing the Ce^i as a power series in et'. 



Na 



a, 



'-'712*^5 ■ 



(14) 



The matching procedure then determines the (in prac- 
tice, much smaller set of) coefficients c„i, rather than 
the Ce^i directly, which eliminates ill conditioning B. 

To extract an A^ x iV coefficient matrix, where A^ — 
Nd + Nc, we need only A^ of the i equations dlC). Alter- 
natively, one may use all i equations as in Poet^ | . In this 



case, the system of equations is overdetermined. Never- 
theless, a solution can be found by the standard method 
of minimizing the sum of the squares of the residuals 
[the differences between the left- and right-hand sides of 
equations (|lO|)]. Previously |^, we found that the least- 
squares method is generally stabler than keeping any sub- 



set of just A^ equations (10) 



Our numerical method is stable and rapidly conver- 
gent. For a given grid spacing h, we established conver- 
gence in propagation distance by performing the match- 
ing every 40 a.u. until convergence was obtained. At 
each matching radius, both the number of discrete chan- 
nels Nd and the number of expansion functions for the 
continuum Nc were varied to obtain convergence. Fi- 
nally, the entire calculation was repeated for a finer grid 
(using one-half the original grid spacing h). 

The biggest advantage of having a general, propagating 
solution is that once the grid spacing is chosen, a "single" 
calculation is all that is needed to establish convergence 
for the remaining numerical parameters. This is because 
the D-matrix, the calculation of which consumes nearly 
all the computational effort, is independent of asymptotic 
boundary conditions. Thus, in a typical calculation, the 
same D-matrix is used for, e.g., Nc — 0, 1, . . . , 9 while 
Nd runs from 1 to 30. This would have required 300 
completely separate calculations (each taking about the 
same time as our "one" calculation) had we solved the 
original global matrix equations (0). 

We have performed complete calculations for electrons 
colliding with atomic hydrogen at impact energies of 54.4 
and 40.8 eV (total energies of 3 and 2 Ryd, respectively). 
In Table I, we present our calculated cross sections for 
e^ + H(ls) -^ e^ + H(ns),?T. < 8. The grid spacing is 
h — 1/5 a.u. (results using one-half this spacing differed 
by less than 0.1% for discrete excitations and 0.5% for 
elastic scattering). One of the advantages of our direct 
approach is that we are able to obtain the amplitudes 
for higher-level (Rydberg) transitions as easily as those 
for low- level excitations, provided the matching radius is 
large enough to enclose the final Rydberg state. This is 
in contrast to some other approaches, such as the CCC, 
which lose accuracy for higher- level transitions. 

TABLE I. Cross sections (vrao; superscripts indicate pow- 



ers of 10) for e- + H(ls) 
of 54.4 and 40.8 eV. 



+ H(ns), for impact energies 



54.4 eV 
Singlet Triplet 



6.47"^ 
4.66"=* 
1.22"^ 
4.92"* 
2.48"* 
1.42-" 
8.89"^ 
5.94-5 



4.07"' 
4.04-^ 
8.39-'' 
3.I3-* 
1.52-* 
8.52-5 
5.27-5 
3.49-5 



40.8 eV 
Singlet Triplet 



8.58-^ 
8.09-^ 
2.15-^ 
8.74-* 
4.41-" 
2.53-" 
1.58-" 
1.06-" 



6.34-" 
5.08-^ 
9.88-* 
3.59-* 
1.71-* 
9.55-5 
5.88-5 
3.88-5 
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FIG. 1. Singlet SDCS (7rao/Ryd) vs. the energy fraction 
thjE for an impact energy of 54.4 eV. The total ionization 
cross sections from FDM, ECS, and CCC are 1.50^^1. 54"^ 
and 1.48^^ ij^a^o), respectively. 



FIG. 3. Singlet SDCS (jrao/Ryd) vs. the energy fraction 
Eft/iS for an impact energy of 40.8 eV. The total ionization 
cross sections from FDM, ECS, and CCC are 1.97"^ 2.04"^ 
and 2.02^'^ (""ao), respectively. 
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FIG. 2. Same as Fig. 1 for the triplet case. The total ion- 
ization cross sections from FDM, ECS, and CCC are 3.11^"^, 
3.39^^, and 3.21~^ ij^a^c^-, respectively. 



FIG. 4. Same as Fig. 3 for the triplet case. The total ion- 
ization cross sections from FDM, ECS, and CCC are 2.47^^, 
2.70"^, and 2.50"^ (ttOo), respectively. 



In Figures 1-4, we present our results (labeled FDM for 
finite-difference method) for the single-differential cross 
section (SDCS). For a total energy of 3 Ryd, 240 a.u. 
proved to be a sufficient matching radius to get conver- 
gence of the SDCS and for E = 2 Ryd, a radius of 360 a.u. 
was required. The SDCS is more sensitive to the number 
of expansion functions for the continuum than the other 
observables, particularly about ei, = E/2. Nevertheless, 
convergence to better than 1% was readily obtained us- 
ing 7-8 functions (the largest discrepancy in the SDCS 
between Nc = 7 and Nc = 8 was smaller than 0.3%; even 
using just 6 expansion function gave results accurate to 
1%). 

Also shown in Figs. 1-4 are the results of convergent 
close-coupling (CCC) calculations [y|. The CCC method 
of Bray g employs a "distinguishable electron" prescrip- 
tion, which produces energy distributions that are not 
symmetric about e^ = E/2. Stelbovics g] has shown that 
a properly symmetrized CCC amplitude yields SDCS 
that are symmetric about E/2 as well as being four 
times larger at £(, = E/2 than those assuming distin- 
guishable electrons. (Note that our singlet FDM results 
at £(, = E/2 are about four times larger than the corre- 
sponding CCC results.) Other than making the energy 
distributions symmetric, it is clear from the figures that 
symmetrization (coherent summation of the CCC ampli- 
tudes Ce^i and Ce^i, where Eq = E — e^, which corre- 
spond to physically indistinguishable processes) will sig- 
nificantly affect only singlet scattering (and then only 
near £{, = E/2), since C^^i is practically zero for e?, > 
E/2. For singlet scattering, the CCC oscillates about the 
true value of SDCS, except near (and beyond) £(, = E/2. 
CCC results for triplet scattering, on the other hand, are 
in very good agreement with our results for < ef, < E/2. 

Some very recent results from Baertschy et al. [|6| have 
also been included in the figures. Baertschy et al. re- 
arrange the Schrodinger equation to solve for the outgo- 
ing scattered wave. They use a two-dimensional grid like 
ours, but scale the coordinates by a complex phase factor 
beyond a certain radius where the tail of the Coulomb po- 
tential is ignored. As a result, the scattered wave decays 
like an ordinary bound state beyond this cut-off radius, 
which makes the asymptotic boundary conditions very 
simple. By computing the outgoing flux directly from 
the scattered wave at several large cut-off radii, and ex- 
trapolating to infinity, they obtain the single-differential 
ionization cross section without having to use Coulomb 
three-body boundary conditions. This method, called 
exterior complex scaling (ECS), has just been extended 
to the full electron- hydrogen ionization problem H. It 
is seen from Figs. 1-4 that the ECS results are in good 
agreement with our FDM results except when the energy 
fraction e^/E approaches or 1. Baertschy et al. M note 
that their method may be unreliable as et approaches 
or E due to "contamination" of the ionization flux by 
contributions from discrete excitations. 



We note also the recent work of Miyashita et al. Q, 
who have presented SDCS for total energies of 4, 2, and 
0.1 Ryd using two different methods. One produces an 
asymmetric energy distribution similar to that of CCC 
while the other gives a symmetric distribution. Both 
contain oscillations. The mean of their symmetric curve 
a.t E = 2 Ryd (40.8 eV impact energy) is in reasonable 
agreement with our calculations. 

In conclusion, we have presented complete, precision 
results for the Temkin-Poet electron-hydrogen scattering 
problem for impact energies of 54.4 and 40.8 eV. It may 
be possible to improve the speed of the present method by 
using a variable-spaced grid, like that used by Botero and 
Shertzer [ |lO[ in their finite-element analysis (this would 
greatly reduce storage requirements as well). Once we 
have optimized our code for this simplified model we will 
proceed to include angular momentum. When angular 
momentum is included, the ionization boundary condi- 
tion is no longer separable and this is the major chal- 
lenge for generalizing the present approach to the full 
electron-hydrogen scattering problem. 
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